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Abstract 

In quantum spin chains, it has been observed that the incommen- 
surability occurs near valence-bond-solid (VBS)-type solvable points, 
and the correlation length becomes shortest at VBS-type points. Be- 
sides, the correlation function decays purely exponentially at VBS- 
type points, in contrast with the two-dimensional (2D) Ornstein-Zernicke 
type behavior in other regions with an excitation gap. We propose a 
mechanism to explain the onset of the incommensurability and the 
shortest correlation length at VBS-like points. This theory can be 
applicable for more general cases. 

Frustration is one of the important problem in many-body physics. Re- 
cently, many things in frustrated systems become clearer with field theoretical 
approaches and numerical calculations (density matrix renormalization group 
(DMRG) etc.). However, relatively little is known about commensurate- 
incommensurate (C-IC) transitions induced by frustration in quantum spin 
systems. Among them, it is observed the C-IC change in the region with 
an excitation gap (therefore, it is not a phase transition) in one-dimensional 
(ID) quantum spin models. 

For example, the C-IC change begins from the Majumdar-Ghosh point 
|]J in the S=l/2 next-nearest-neighbor (NNN) chain 

H = XX^+i + aSjS j+2 ), (1) 



1 



(«d = 1/2 is the solvable point), or from the valence-bond-solid (VBS) point 
in the S=l bilinear-biquadratic chain, 

H = J2(S J S J+1 + a(S J S J+1 ) 2 ), (2) 
j 

{old = 1/3 is the solvable point). 

With a numerical diagonalization, Tonegawa and Harada observed in the 
S=l/2 NNN chain || that the incommensurability of the spin correlation 
begins near the solvable point. However, from the DMRG calculation, Bur- 
sill et al. stated that in the static structure factor (Fourier transform of 
spin-correlation function), the C-IC change begins slightly different from 
the solvable point for the S=l bilinear-biquadratic chain [Q, and for the 
S=l/2 NNN chain ||. In contrast, using the DMRG, Schollwock, Jolicoeur, 
and Garel || have closely investigated the correlation function of the S=l 
bilinear-biquadratic chain (0), not the structure factor, and they found that 
the incommensurability begins right at the VBS point. Adding this, they 
have pointed out the singular behavior of the modulation wavenumber. 

Secondly, it has been noticed that the correlation length is shortest at the 
VBS-type point. 

Thirdly, at the VBS-type point, the spin-correlation function decays purely 
exponentially, in contrast with the 2D (1+1D) Ornstein-Zernicke type (that 
is, 2D Green function for the free massive particle) behavior in the other 
region with a gap. 

Although Schollwock et al. have tried to explain these features on the 
analogy of the ID classical frustrated Ising spin with finite temperature, 
several features are different from the quantum spin chain case, which we 
intend to explain in this study. After that, Kolezhuk, Roth, and Schollwock 
applied their idea to the S=l NNN chain 0, [J, where they called the shortest 
correlation length point as the "disordered point", generalizing the VBS-type 
point with the above features. 

Recently, Fath and Siito proposed an effective field theory Q to explain 
these features in the C-IC change of the S=l bilinear-biquadratic chain. Al- 
though their results seem natural, they have not derived their effective field 
theory directly from the lattice Hamiltonian. And from their theory, the C- 
IC change near the Majumdar-Gosh point cannot be explained. Therefore, 
it is meaningful to construct a theory on firmer ground. Considering the 
analyticity of the structure factor S(q) around the disordered point, where 



2 



several exact results are known, we propose another approach to this prob- 
lem. Since, what we observe experimentally is not Hamiltonian itself, but 
structure factors or other quantities, thus our approach may be more natural 
in a physical sense. 

We summarize the results by Schollwock et al. as (see also Fig.jl]) 




Figure 1: Singular behaviors of correlation length (a) and wave number (b). 

1. In the commensurate region (a < old), the correlation function behaves 
as the 2D Ornstein-Zernicke-type 

(S S n ) « (-lrM-^expHnl/^a)), (3) 

where far enough from the disordered point old- 

2. In the incommensurate region (a > old) , 

(S S n ) « (-l) n cos(g(a)n)|n|- 1 / 2 exp(-|n|/e(a)), (4) 
where far enough from the disordered point old- 

3. At the disordered point old, 

(S S n ) w (-l) n exp(-\n\/^a D )). (5) 



4. From the commensurate side, 



da 

5. From the incommensurate side, 

dq(a 



\OLD 



-oo. (6) 



(7) 



da 

6. From the incommensurate side, by fitting as 

q(a) oc (a — ao) ff , (8) 

we obtain a ~ 1/2. 

7. Neither the exp(— |n|/£) nor the |n| -1 / 2 exp(— |n|/£) form cannot de- 
scribe the correlation function near the disordered point op. 

Note that, strictly speaking, the 2D Ornstein-Zernicke form is the modified 
Bessel function |ID[ , 



G(fl) = / expwy = 2lKMnm (9) 



which behaves |n| 1//2 exp(— |n|/£) in large distance. And static correlation 
function is 

C (0 1 .,)=J^M i(}1 . (10) 

To explain these features, we consider the static structure factor in the 
complex q plain, 

s( q ) = ^Mkn)(-i) n (s s n ). (ii) 

n 

The simple minded Fourier transforms of the correlation functions are 
1. At the disordered point, 

S(q) oc — L-j. (12) 
q l + rrr 
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2. In the commensurate region, 



S(q) oc 



1 



(13) 



\/q 2 + m 2 



3. In the incommensurate region, 



S(q) oc 



1 



1 



(14) 



+ 



(q — q(a)) 2 + m 2 



(q + q(a)) 2 + m 2 



However, with the above forms, one cannot connect these three region. In 
fact, in terms of the singularity in the complex q plain, there are poles at the 
disordered point, in contrast with the branch-cuts in the other (commensu- 
rate, incommensurate) regions. 

In order to unify these three region, we reconsider the relation between 
the pole and the branch-cut. A pole and a branch-cut can be deformed each 
other as following, 



d — 0, and two branch-points at z = ±iJ\d\ for d < 0. Note that f(z) is an 
odd function f(—z) = —f(z), considering the analytic continuation around 
the branch-cut. 

Besides, to relate it with physical phenomena, we consider the following 
assumptions: 

1. S(q) is analytic except several algebraic singular points. 

2. S(q) is real on the real axis (S(q) = S(q)). 

3. From parity, S(q) = S(—q). 

4. S(q) is varied continuously with a parameter a. 

5. At the disordered point ap, S(q) is described with two poles. 

6. Singular points and branch-cuts should not cross the real axis. 

From these assumptions, we can obtain the simplest singularity structure 
as Fig. H, that is, there are two poles symmetric with real axis, or four 
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branch-points continuously deformable from them. The static structure fac- 
tor, originated from the singular terms, is given as 

S S ing(q) = Af(q + mi)f(q - mi) 

= A[q A + 2(m 2 -d)q 2 + m 4 + d 2 + 2m 2 d]~ 1 / 2 , (16) 

where real parameters A, rh, d may depend on a, especially d oc a — a D . 
Although this result for the static structure is the same as in M, our theory 
is based on the analyticity of S(q) and the exact S(q) at the VBS point. 
Note that without the assumption for algebraic singularities, there remain 
other possibilities 

i 

S sing (q) = A—[f(q + mi)-f(q-mi)] 

: cos ( - arctan( — ^ — -) \}7) 



™J(q 2 -rh 2 -d) 2 + 4q 2 m 2 \ 2 q 2 - m 2 - d 



S sing (q) = A^[f(q + mi)+f(q-mi)]. (18) 
The assumption for algebraic singularities physically relates to the locality of 
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Although it seems reasonable four branch-points in the incommensurate 
region, how do we interpret four branch-points in the commensurate region? 
In the large distance, or the small wave number q, to S(q) on the real axis 
mainly contribute the two branch-points closer to the real axis, thus it seems 
only two branch-points when < 1 — \d\/m 2 <C 1 . 

When we set d « c\{a — an)] C\ > in eq. (|T6"D, we can explain system- 
atically the singular behaviors of £(a),g(a) near the disordered point «». 
And we can calculate the analytic form for S(q) in the neighborhood of «o. 
In addition, near a D , m may behave asm~ m + c 2 (a — a_o). It should be 
C2 < 0, since the correlation length increases in the incommensurate region. 

Generally, the requirement for the amplitude is A ^ 0, since for A = 
the correlation function becomes perfectly zero. However, at the Majumdar- 
Ghosh point in the S=l/2 NNN chain, the correlation function behaves as 

(S S n ) = 0, (for\n\ > 2), (19) 

which can be interpreted as the A(a) oc (a—ap) case or the A(a) oc (a— o>d) 2 
case. The latter possibility comes from that spin variables may depend as 
S n oc (a — ct£)). 
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Figure 2: Continuous deformation of branch-cuts and poles. They correspond 
to commensurate region (a), disordered point (b), incommensurate region (c), 



Until now, we have neglected the lattice structure and the terms which 
are regular in the whole complex q plain. To include the lattice structure, 
we should require S(q) = S(q + 27r). Thus, the full static structure factor is 

oo 

S(q) = Yl S sing {q + 2nj)+S reg {q), (20) 

j=-oo 

where 

oo 

Sreg(q) = J2 a jCOs(jq). 

j=-oo 

Note that to calculate the spin-correlation function, it is enough to consider 
the singularities of S(q) within the Brillioun zone, because of periodicity. 

In summary, the C-IC change can be regarded as a fusion-fission of singu- 
larities in the complex q plain. The Lifshitz point, where the static structure 
factor changes from a single peak to double peaks, is a natural consequence 
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of this fusion- fission of singularities || 0, |3| (see in ( |T6"D , 

g^M ( g = 0) ex A (m 2 - d) (21) 

)• 

The argument in this paper should be checked in the DMRG calculation 
near the disordered point. After checking the expectation for the static struc- 
ture factor, one can extend it to the dynamical structure factor. Since the 
VBS-type states || are obtained not only in the ID quantum system, but 
also in the 2D and 3D quantum systems, it is interesting whether such C-IC 
changes may occur at the VBS-type points in higher dimensions. 
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